MULTIVARIATE POLYNOMIAL INTERPOLATION AND SAMPLING IN 

PALEY- WIENER SPACES 



B. A. BAILEY 

Abstract. In this paper, an equivalence between existence of particular exponential Riesz bases 
for multivariate bandlimited functions and existence of certain polynomial interpolants for these 
bandlimited functions is given. For certain classes of unequally spaced data nodes and corresponding 
£2 data, the existence of these polynomial interpolants allows for a simple recovery formula for 
multivariate bandlimited functions which demonstrates L2 and uniform convergence on R d . A 
simpler computational version of this recovery formula is also given, at the cost of replacing L2 
and uniform convergence on R d with L2 and uniform convergence on increasingly large subsets of 
R d . As a special case, the polynomial interpolants of given £2 data converge in the same fashion to 
the multivariate bandlimited interpolant of that same data. Concrete examples of pertinant Riesz 
bases and unequally spaced data nodes are also given. 



1. Introduction 

Approximation of bandlimited functions as limits of polynomials has a long history, as the following 
question illustrates: if (sinc7r(- — i n ))nez forms a Riesz basis for PWi—^^u what are the canonical 
product expansions of the biorthogonal functions for this Riesz basis? The first results along these 
lines were given by Paley and Wiener in [8], and improved upon by Levinson in [5, pages 47-67]), 
while Levin extends these results to different classes of Riesz bases in pQ. A complete solution 
is given by Lyubarskii and Seip in [6] and Pavlov in |9]. In particular, they prove the following 
theorem: 

Theorem 1.1. Let (t n ) n C R, where t n 7^ when n ^ 0, be a sequence such that the family of 
functions (sinc7r(- — t n )) n is a Riesz basis for PW_ T j, then the function 

s W = ^ (,-*,) n (!-£) 

{tn ■ \t n \<r,n^0} 

is entire, where convergence is uniform on compacta, and the biorthogonal functions (G n ) n of 

(sinc7r((-) - t n )) n 

are given by 

G n (z) - S{Z) 



(z-t n )S'(t n ) 

The following is a readily proven corollary of Theorem 11.11 
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Corollary 1.2. Let (t n ) n C R and {Gk)k be defined as in Theorem \l.ll then for each k, there 
exists a sequence of polynomials ($jv fc)jv such that 

1) ®N,k{tn) = G k (t n ) when \t n \ < N. 

2) limjv->oo &N,k = Gk uniformly on compacta. 
Corollary 11.21 raises two questions: 

1) Does every multivariate bandlimited function, (not just biorthogonal functions associated with 
a particular exponential Riesz basis), have a corresponding sequence of polynomial interpolants? 

2) If such polynomial interpolants for a multivariate bandlimited function exist, can these inter- 
polants be used to be approximate the function in some simple and straightforward way? 

Let (in) rae z d C M rf be a sequence such that the family of exponentials ( e ) ne gd i s a uniformly 
invertible Riesz basis for L2([— vr, Tr] d ) (defined in section 4). Under this condition, Theorem 14.21 
answers the first question affirmatively by showing that multivariate bandlimited functions can be 
approximated globally, both uniformly and in L2, by a rational function times a multivariate sine 
function. Stated informally, 

T / s SINC(7rt) 
Qd,e{t) 

where (^JJfeN is a particular sequence of interpolating polynomials and (Qd,i)d,t is a sequence 
of polynomials which removes the zeros of the SINC function. This gives a partial answer to the 
second question, but the fraction in the approximants above becomes more complex as £ increases. 
Theorem 15.11 gives a more satisfactory answer to question 2) by using 

_ x^N 1 II* II if 

e Z^fc=l k {2k-\) (£+i/ 2 )2fc-l ^ £>0 

in lieu of the SINC function in expression (fTJ. The exponent in the above expression is simply a 
rational function of I. This simplicity necessitates replacing global L2 and uniform convergence 
with a more local (though not totally local) convergence. Corollary 15.21 is of particular interest as 
a multivariate analogue of Theorem II. 1| stately informally as 

f(t) ~ ^(t), £ > 0. 

The author is unaware of any other multivariate polynomial approximation theorem which applies 
to exponential Riesz bases which are not necessarily tensor products of single- variable Riesz bases, 
or that demonstrate convergence stronger than uniform convergence on compacta. As a note, The- 
orems 14.21 15.11 and Corollary 15.21 do not, at this point, recover Corollary 11.21 in its generality of 
allowable sequences {tn)n C R; however, the comments above show that their value is primarily 
due to their multidimensional nature and convergence properties. 



This paper is outlined as follows. Section 2 covers the necessary preliminary and background 
material regarding bandlimited functions, and section 3 outlines some basic properties of uniformly 
invertible operators. Theorems 14.21 and Theorem 15.11 are proven in sections 4 and 5 respectively, 
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along with pertinant corollaries. Section 6 gives explicit examples of sequences (t n ) n£ %d to which 
Theorems 14.21 and 15.11 apply. Section 7 (as an appendix) addresses the optimality of growth rates 
appearing in Theorem 15.11 

2. PRELIMINARIES 

Definition 2.1. A reproducing kernel Hilbert space is a Hilbert space H of functions on X such 
that there exists K : X x X — > C satisfying the following: 

1) For all y G X, K(-,y) G H. 

2) f(x) = (f(-),K(-,x)) for all x G X and / G H. 

Definition 2.2. A Riesz basis for a Hilbert space H is a sequence (/ n )neN which is isomorphically 
equivalent to an orthonormal basis of H. Equivalently, a Riesz basis is an unconditional Schauder 
basis. 

If (/n)neN is a Schauder (Riesz) basis for a Hilbert space H, then there exists a unique set of 
functions (/^) n eN ; (the biorthogonals of (f n )n<=N) such that {fn, fm) = ^nm- The biorthogonals 
also form a Schauder (Riesz) basis for H. Note that biorthogonality is preserved under a unitary 
transformation. 

We use the (i-dimensional L2 isometric Fourier transform 

*(/>(■ ) = P-V. AO."" -M. / € L 2 (R d ), 

where the inverse transform is given by 

^'(/)(- ) = P.V. -L^ I ! 6 L 2 (M«). 

Definition 2.3. We define PW[_ W)7r] d := {/ G L 2 (R d )|supp(.F- 1 (/)) C [-vr,7r] d }, with the inher- 
ited L2(M d ) norm. Functions in PWy_ n ^] d are a l so called bandlimited functions. 

Here are facts concerning PWi w ^w which will be used ubiquitously. 

1) PWi^jd is isometric to L2([— n, 7r] d ) by way of the Fourier transform. 

2) PWr^id consists of entire functions, though in this paper we restrict the domain to M. d . 

3) PWr^id is a reproducing kernel Hilbert space with reproducing kernel 

K(x,y) = SINCvr(x-y) 

where 

SINC(x) := sinc(x(l)) • . . . • smc(x(d)), sinc(x) := 

x 

4) (SINC7r((-) — ")) neZ d i s an orthonormal basis for PWr^^u- This follows from 

^(^e ir (-))(t) = sincTrfi - r). 
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5) In PWt^^id, L 2 convergence implies uniform convergence. 

6) If / € PW { _^ ]d , then 

lim f(x) = 0. 

This follows from the d-dimensional Riemann-Lebesgue Lemma. 

7) The following result [10} Theorem 19.3] due to Paley and Wiener characterizes single-variable 
bandlimited functions. 

Theorem 2.4. A function f is in PW r |_ I)I j if and only if the following statements hold. 

1) f is entire. 

2) There exists M>0 such that \f(z)\ < Me^ for zeC. 

3) f\ R EL 2 (R). 

3. Uniform invertibility of operators and Riesz bases 

Definition 3.1. Let A : ^(N) —> ^(N) be an onto isomorphism. Regard A as a matrix map with 
respect to the unit vector basis of ^(N). Let be the orthogonal projection onto the span of the 
first k terms of the unit vector basis. If 

(2) sup||(vr fc .^TTfc/r 1 !! < oo, 

for an increasing sequence (kj)j^, then A is said to be uniformly invertible as a matrix map with 
respect to the projections ("Kk^jeN- The terms in inequality ([2]) should be interpreted as standard 
matrix norms and inverses of finite dimensional matrices. 

Let S be an orthonormal basis for a Hilbert space H . Let (SVi)neN be a sequence of sets such that 

1) + Si C S 2 C • • • C S, and 

2) IT=1 Sn = S. 

Define Pi to be the orthogonal projection onto sp&n{ep.}e h ESc Note that 

(3) lim P^x = x, x £ H. 

Linearly order S = (e n ) ng N such that, if e n G \ S^-i, and e m £ S&_i, then m < n. 



Definition 3.2. Let (vk)k£N be the unit vector basis for ^(N) and define <f> by 4>ek = Vk- Let 
L : H — > H be an onto isomorphism. L is said to be uniformly invertible with respect to the 
projections (Pi)i£fq if (frLcj) -1 is uniformly invertible as a matrix map with respect to the projections 

(^l)- 

We define the following notation: 

(4) (P,LP,)- 1 :=(vr |s£| (0L0- 1 )^ | )- 1 . 

By saying PgLPi is invertible, we mean that the right hand side of equation (jH is well defined. If 
L is defined on span(e n ) n6 N, but perhaps not on H, we define tl PiLPi is invertible" in the same way. 
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If L is an operator on H (perhaps densely defined), and (Pe)eeN is a sequence of projections denned 
above, define the operator L# = LPe + I — P#. 

Definition 3.3. Let (vk)keN be a Riesz basis for H. We define (ufc)fceN to be uniformly invertible 
with respect to the projections (Pe)e<^n if the corresponding isomorphism Le^ = vj. is uniformly 
invertible with respect to the projections {Pe)e&i- 

We can now state and prove the following lemmas: 

Lemma 3.4. Let (e n ) n£ N be an orthonormal basis for H , let (/ n )neN C H, and let Pi be the or- 
thogonal projection onto span(e n ) n <^. Define L : span{e n } n£ N — > H by Le n = f n . For each £ > 0, 
the following statements are equivalent: 

1) {fn)n<t U (e n ) n>i is a Riesz basis for H. 

2) P(LPi is invertible. 

Proof of Lemma \3.4\ 1) ==^ 2): From the definition of hi we know that it is an onto isomorphism. 
This yields P e = P e LP e Lj\ implying P e = (PtLP £ )(P e LfP e ). 

2) ==> 1): Let Ae be the unique square matrix such that PeLPgAe = AePeLPe = Pe. We need to 
show that Li is an onto isomorphism. 

First we show that Le is one to one. Say = L^x = LP^x + {I — Pg)x, then = PeLPex, so that 
= AiPpLPix = P^x. We conclude that x = (I—P^jx. This implies = Le(L—Pe)x = (J— Pg)x = x. 

Next we show that Le is onto. Note Le(L — Pg)x = (I — Pg)x, so we only need to show that for all 
x, Pix is in the range of Le- Define 

y = P e A e P f x + P e x - LP e A f P e x. 

We have Le(PeAePex) = LP^A^Ppx and 

(5) L e (P e x - LP e A e P e x) = LP t x - L(P e LP e )(A e P e )x + P e x - LP e A e P e x 

= P lX - LPeAePex, 

from which Ley = P^x. 

Le is a continuous bijection between Hilbert spaces, and hence is an onto isomorphism by the open 
mapping theorem. □ 

Lemma 3.5. Define L as in Lemma \3J\ For all I > 0, Le is an onto isomorphism iff it is one to 
one. 

Proof of Lemma \3.5[ We only need to show that PfLPi is one to one on PeH when Le is one to 
one on H, and apply Lemma 13.41 Let (PeLPe)Pex = 0. We have 

L e [Pex - (L - P e )LP e x)] = LePex - Le(L - Pe)LPex 

= LePex - (I - Pe)LPex 
= LePex — LPex = 0, 
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where the last equality follows from LpPp = LP?. Since Lp is one to one, we have that Pgx = 
(I - Pi)LPiX, so that PiX = 0. □ 

Lemma 3.6. Let (e n ) ne N be an orthonormal basis for H , and (/ n )neN be a Riesz basis for H, and 
(^)feN C N be an increasing sequence. Let Pi be the orthogonal projection onto span{e n } n <fc, ; then 
the following are equivalent: 

1) The operator L is uniformly invertible with respect to (P^)^gN- 

2) For all £ > ; Lp is an onto isomorphism, and sup£ gN HL^H ^ 00 • 

Proof of Lemma \3.b\ 1) ==> 2): By Lemma \3A\ we only need to show that that 

sup || LJ" || < oo. 

This follows from the identity 

(6) L7 1 = [/-(/- Pp)L\{PpLPp)- 1 + L - P t , 
which is hereby demonstrated: 

(7) [J — (J — Pp)L]{PpLPp)~ 1 + L — Pe = [I-tf-PfiLlPtL^Pt + I-Pt 

= PeL e x Pp — LPpL e 1 Pp + / 
= {L -L)PpLj 1 Pp + L. 

We have (I — L)Pg = L — Lp, so 

(8) [I - (I - P^PeLPi)- 1 + I - Pi = (L-L)Lj 1 P i + L 

= LfPe-Pe + I. 

Noting that Lp(L — Pp) = L — Pg, we obtain L^ l Pg — Pp + L = L7 1 , which proves the identity. 
2) => 1): Noting that (PpLPp)' 1 = PpL~ x Pp yields the result. 

□ 

4. The first main result 
We begin with some necessary definitions: 

Definition 4.1. Define C i>d = {-£,..., £} d , and e n (x) = j^y^e i{ - x ^ for n £ Z d . Let Pp : 
L-2{[— 7r ] d ) L/2{[— 7r, ir] d ) be the orthogonal projection from L2Q— n, vr] d ) onto span(e n ) ng c f d - 

Let (/n)nez d De an exponential Riesz basis. In the following sections, we abbreviate the statement 
"(fn)nez d 1S a uniformly invertible Riesz basis with respect to the projections (PejeeN defined in 
definition 14. IT by "(/n) n ez d is a uniformly invertible Riesz basis". 

To avoid confusion of indices, we write t £ M. d as t = (i(l), • • • , t{d)). 
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For £, d G N define the multivariate polynomial 



Here is the first main result of this paper. 

Theorem 4.2. Let (t n ) n€Z d C R d , and define f n {x) = ,» \ d / 2 ^ x,tn ^ for n G Z d . Lei {f n ) n ez d ^ e 
a uniformly invertible Riesz basis, then for all f E PWi^^id, there exists a unique sequence of 
polynomials (\E^)^eNj ^ : M. d — > R, such that 

(a) has coordinate degree at most 2£. 

(b) ^ e (tn) = f(t n ) for all n G C tjd . 

(c) f(t) = lim^oo ^e(t) S q^^(^ > where the limit is both L2 and uniform. 

Note: The expression in statement (c) of Theorem 14.21 has removable singularities, but these can 
be evaluated by 

.. sinc(vri) (£\) 2 

teq^- « + n)i(«-»)i - " e { -'-- e} - 

The proof of Theorem l4.2l requires several lemmas, beginning with the following equivalence between 
the existence of particular Riesz bases and a polynomial interpolation condition: 

Lemma 4.3. Let {t n ) n&d C R d where f n (x) = j^m^ x ^ and e n (x) = (^372 e^' n > . The 

sequence (fe >n )nez d '■= (fn)nec e , d U (en)n<£c e . d ™ « R^sz basis for L 2 ([-tt, 7r] d ) #f £/ie following con- 
ditions hold. 

1) For allneC £jd , t n <£ (Z\{-£,--- ,£}) d . 

2) For any sequence (c k ) k& c ld , there exists a unique polynomial with coordinate degree at most 
2£ such that \E^(ifc) = c k for k G C^. 

Proof of Lemma \4-3[ Suppose that the sequence (fn) n <=z d * s a Riesz basis for L2([— vr, ir] d ). We 
compute the biorthogonal functions of {fe,n) n GZ d when n G Ci& 

fin = ^2(ft^ e k)e k = ^2 {fl n ,e k )e k + ^2 (fl n ,e k )e k 

k& d keC tid HCi,d 

keC e>d 

Passing to the Fourier transform and defining G^ n = J~(f% ), we have 



(9) G e>n (t) = Yl G ttn (k)SmCn(t-k) 



keC ed 



sr c,,„(t)(-i)*<»+ ■■ ■whm t{d) \ 
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Denote the k th summand in equation Q by A^, then 

n\=-* (/(i) .y,)-...-ir, . m-u) 

Hn ' k ~ ^ k w\- t (*d)-ii)-...-n^ (*(co-i-) 

ji^fc(i) i d #fc(d) 

G < , n (fc)(-i)*( i )+-+*( d )t(i) • . . . • nV-* (*(i) ./!)••••• nV-« - ^) 
nl=-^(i)-ii)-----nl=-^)-id) 

G<,n(fc)^(-l)^ + '" + ^nn=-« (*(1) ./:)••••• n^=-^ (t(«0 - 3d) 
giffcgj JdffcW 

nU(>-f)--nU(i-f) 

_ Pl,n,k{t) 

where p£ ini fc is some polynomial with coordinate degree at most 11. Substituting into equation 
we obtain 

Gi,n(t) = ( ^2 Pi,n,k{t)\ -77-7^- : = ^,n(*)" SE ^ ' ' 



fceQ >d 



Qd,£(t) * v ' Qd,*(t) 



where 4>£,n is a polynomial having coordinate degree at most 11. The fact that each zero of sinc(7rz) 
has multiplicity one implies that the zero set of S Q d °ff (which is entire) is (Z \ • • • , l}) d C C. 
Using that SINC-7r(x — y) is the reproducing kernel for PW^^d, we see that 

for all n,m £ ^. This yields that 



This shows 4>e jn (t n ) 7^ and / for n G and that i n ^ (Z\ {— £, ■ ■ ■ ,£}) d (statement 

1) of Lemma 14.31). 



For n, m £ Cg^, n 7^ m, 

n s~i (4- \ J. I x/SINC(7Tt)\ 

— (j£,n{tm) — <Pe,n{tm) I — ^ 777— I 

V Qd,£\t) ' 

We conclude that 

I Qd,i(tn) 

4>£,n(t m ) = 



SINCvri n ^ 0' n — m 

0, n^m 



for n,m £ Cg d From this, the "existence" part of statement 2) in Lemma 14.31 readily follows. Re- 
stated, the evaluation map taking the space of all polynomials of coordinate degree at most 11 to 
jj(2£+l) j g on to. These spaces have the same dimension, hence the evaluation map is a bijection, 
which completes the proof of statement 2). 
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Suppose that 1) and 2) hold. For n G Cp^-, let p^ n be the unique polynomial of coordinate degree 
at most 2£ such that pi, n (t m ) = 5 nm for m G C^. Define 

(10) *'' n(t) = Q rf/ (t)SINC^/ ^ (t) - 

Partial fraction decomposition can be used to show that G PW^^^w. For n, m G C^d, we 
therefore have 

< m = (^,„(-),SINC7r((-) -i m )> = (J^-H^Jm). 
Let be defined as before. Let / = X^nez d c n e n such that Lp(f) = 0, then 

— ^ ^ C-nfn ~t~ ^ ^ C n 6 n . 

If, for each n G Ci : d we integrate the above equation against J r_1 (<I>^ in ), we see that Ck = for 
& ^ C^,di so that Cfc = for all G is one to one, so by Lemma [3.5l it is an onto isomorphism 

from L 2 ([-iT,7r} d ) to itself. □ 

Proof of (a) and (b ) of Theorem \4-2\ Lemmas 13.61 and 14.31 imply the existence of a unique sequence 
of polynomials satisfying statements (a) and (b) of Theorem 14. 2\ namely, 

Yl /(^Kn(t), 

where pp tn is defined as in the proof of Lemma 14.31 



□ 

It remains to show that this sequence of polynomials satisfies the statement (c) of Theorem 14.21 
Proposition 4.4. The following are equivalent. 

1) L is uniformly invertible with respect to the projections (Pe)eeN- 

2) For all x G H, lim^ 00 (L|)- 1 (I - P £ )x = 0. 



Proof of Proposition \4-4\ It is clear that 1) implies 2). For the other direction, note that the 
equality L* e = PpL* + I — Pp implies that 

(11) /=(L|)- 1 P,L* + (L|)- 1 (/-P,). 

This implies that (L|) _1 P^ is pointwise bounded. Together with the assumption in 2), this implies 
is pointwise bounded, hence norm bounded by the uniform boundedness principle. This 
yields uniform invertibility of L. □ 

Lemma 4.5. The following are equivalent: 

1) For all g G -^Q - tt, 7T ] d ), we have 

(12) g= \\m{L})- l P,L*g. 

£—>oo 

2) L is uniformly invertible. 

Proof of Lemma \4-5[ Recall equation (llip and apply Proposition 14.41 □ 
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Proposition 4.6. Statement (c) of Theorem\4-^\ is true iff 



(13) 



0= lim \f(tn) 



neCtd 



SINCvrtn 



QdAtn) 



lim Si >d , f G PWr id. 



Proof of Proposition \4-6\ Note that Le n = f n implies that /* = (L*) l e n . Similarly, f^ r 
(L})- l e n . Given / G PW^^a, let g = F^ 1 (f). Equation JT2D shows: 



p-\f) = lim (L,*)- 1 V (L*g,e n )e n = jimmy 1 V (gj n )e r 



lim V (g,f n )fl n = lim V /(i„)/| n - 



Passing to the Fourier transform, we have 



(14) 



/= lim J2 f^Hftnl f£PW { ^ ]d , 



where convergence in both L2 and uniform. The values of a function in PW^_ n on the set 
(t n ) n£ id uniquely determine the function. This and equation (fl~0|) show that 

Q d/ (t n )$WCTTt 



Hfln){t)=G l , n {t) 



Qd^SWClTtn 1 



e,n(t), n G Ci t d- 



This implies that 



SINC(TTt) 

Qd,i(t) 



„Z7^ y Vd,nr) 

E /(*»)^t^»)(*)- 



QdAtr 



Combined with equation (|14p . we see that statement (c) of Theorem 14.21 holds iff 

SINC^t n " 



0= lim /(*») 



neCp 



1 



QdAtn) 



P(ft n ), fePW,^d, 



where the limit is in the L2 sense. Passing to the inverse Fourier transform, statement (c) holds iff 

SINCvrt, 



(15) 



0= lim (LI)- 1 /(* 



1 



QdAU 



In addition to having uniformly norm bounded inverses, Lg is pointwise bounded, so there exists 
m,M > such that for all £ G N, / G PIT^, m||/|| < || (L*)- 1 /)! < M||/||. This, combined 
with equation (fTo]) , proves the proposition. □ 

Proposition 4.7. TTie following are true: 



V SUPa.gR sup ter 



sinc(7ra;) 



1. 
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dl -k- < ^3} forte N, then 



2) Define A e4 = jn € 
(16) 



< 1 - S ' NC <""> < 1 

Qd,e(t n ) 



-d(l+2) 

e * 4/3 -i , neAy. 



Proof of Proposition \4~/7\ For 1), the identity 



sinc(vrt) = (l 



k=l 



k 2 



implies 
(17) 



sinc(7ri) 

QiAt) 



n (1 



k=i+l 



k 2 / 1 



where convergence is uniform on compact subsets of C. Fix £ E N. If t G [0,^+1], then 
Note that |Qi,£(i)| = Y\k=i ( vi ~ 1 ) * s increasing on (£ + 1, 00). If t E {£ + 1, 00), then 



sinc(-7rt) 
Qi,£(t) 



< 1. 



sinc(-7rt) sin(7rt) 



< 



1 



Computation yields 



{2£+l)\ 



£l(£ + l) 



so 



sinc(7rt) 



(£\) 2 

< ,\ < 1. 



Q M (t) tt(2^ + 1)! 



Observing that '^jnr is even proves 1). 



For 2), let t G R such that 
(18) - log 



£+1 
sinc(7rt) 



< ^3,then0< and 



fc=£+l 
00 00 1 



fe=f+l j = l J 



Basic calculus shows that 



I27 



< 



j=l J Jfe=£fl 



+ 



fc=£+l 



fc2j ^ (£+l)2j ( 2 j -!)(£+ l)2i-l" 



Equality ([To]) implies 



(19) 



log 



/sinc(7ri)\ 1 / t \ 2 J .„ 1 ft 
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If n € A Ad , then for each 1 < k < d, > jtjg^, so that 

/ SINC(7rt w ) x _ A / sinc(7rt w (fc)) \ + 2) 

° g V Q^(i n ) / ^ ° g V Q M (t„(fc)) / £ 4 / 3 - 1 • 

Statement 2) of Proposition 14.71 follows readily. □ 

Proof of statement (c) in Theorem \4-S\ Proposition 14.71 gives the following: 



st* < ( E + E 

neA £ , d n& d \A i>d 



ntf 



SINCvri r 1 2 



Qd,i{tn, 

/ -djt+2) \ ^ 

(20) < h-e?^rrj £|/(t n )| 2 + ^ 4|/(t n )| 2 . 

n ^ d n: ^<||t»||<» 

Now (/(^n))nez d £ ^2(2^) implies that lim^oo S^d = 0, so that by Proposition 14. 6( statement (c) 
in Theorem 14.21 is true. □ 



5. The second main result 
Theorem 14.21 can be simplified. The function 

SINC(TTt) 

becomes more computationally complex for large values of I. If, at the cost of global L2 and uniform 
convergence, we adopt an approximation 



C-ik=l k(2k-l) 



2 k 
■2k 



(21) SINC(Trt) ~ Q d ,i(t)e ^ fe=1 k{2k - 1] C<+i/a) 

we bypass this difficulty as the exponent of the above quantity is simply a rational function of 
1 > 0. This is precisely quantified in Theorem 15.11 which is the second main result of this paper. 

Theorem 5.1. Let (t n ) Z d C M. d be a sequence such that the associated exponential functions are 
a uniformly invertible Riesz basis for L2Q— 7r, 7r] d ). For N G {0,1,2,...}, A > 0, let Egjf^ = 

2N+1 2JV+1 

[—A(£ + l/2)a^+a,A(^ + l/2)aJV+a]. Let f G PWr^^w where {^e)e is the sequence of interpolating 
polynomials from Theorem \4-ty and let N > 0. Let f G PWr^^w, and define 

I f e (t) = V t {t)e k=1 k(2k - 1] (WW- 1 , 

then 



(22) Mm f{t)-I f/ {t) 



= 0, 

L 2 <,{E ltN>A y) 
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and 
(23) 



lim 

l— >oo 



f(t)-If,i(t) 



£oo((£f,iV,A) d ) 



0. 



If N = in Theorem 15, 1\ we have the following extension of Corollary 11.21 to arbitrary multivariate 
bandlimited functions (at the expense of introducing uniform invertibility) : 



Corollary 5.2. For all f G PWr^^id, we have 



(24) 

and 
(25) 



lim 

£— >oo 



lim 

t— T oo 



f{t) - Mt) 



f(t) - Mt) 



L 2 ([-A(£+l/2) 1 /2 iJ 4(£+l/2) 1 /2]d) 



L 00 ([-A(^+l/2) 1 /2,A(^+l/2)V2]d) 



0. 



0. 



It is evident that if {t n ) n &i d C M rf is any subset such that the associated sequence of exponen- 
tials (f n ) n& z d ^ s a Riesz basis for L2([— it, vr] d ), then the map / h-> (/(in)) ne z d ^ s a bijection from 
-f > W A [- 7r ,7r] d to ^2{^ d )- This allows for a nice interpretation of Corollary 15.21 Given a sequence 
(t n ) ngZ d C R rf (subject to the hypotheses of Theorem 15. ip . and sampled data ((i n , c n )) ngZd where 
(c n ) ngZ d G £2(% d ), then a unique sequence of Lagrangian polynomial interpolants exists, and in the 
appropriate limit, converges to the unique band- limited interpolant of the same data. 

When N = 1, we have a sampling theorem with a Gaussian multiplier: 

|| t ||2 

Compare Theorem 15.11 with Theorem 2.6 in [2], which is a multivariate sampling theorem with a 
Gaussian multipler with global L2 and uniform convergence. Also compare Theorem 15.11 with The- 
orem 2.1 in [11] . which, when d = 1 and the data nodes are equally spaced, gives another recovery 
formula involving a Gaussian mulitplier in the context of over-sampling. 

The proof of Theorem 15.11 relies on two lemmas: 

Lemma 5.3. Let d G N ; N G {0, 1, 2, . . .}, and A > 0. There exists M > such that for sufficiently 
large t, and any t G {E^ n A) d , we have 



"Z^fc=i 



< M(£ + l/2)"^T|SINC(7rt)|. 
Lemma 5.4. For all f G PWr^^id and N G {0, 1,2,...}, we have 



ntii 2(Ar+1) 

11 "2(iV+l) 

g (£+l/2)2JV + l(JV+l)(2JV+l) SINC (7Tt ) 



lim sup 



|| t ,|2(iV+l) 

11 H 2(iV+l) 

g (£+l/2) 2JV + 1 (JV+l)(2JV+l) 



1 /(*) 



The proof of Lemma 15.31 relies on the following proposition. 
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Proposition 5.5. If f : (0, oo) — > (0, oo) is convex, decreasing, differentiable, and integrable away 
from ; then 



1 °° roo 

-f'{l+l/2)< V /(*)- / f(x)dx<0. 
4 ,/T;-, Je+i/2 



k=i+i Je + 1 / 2 

Proof of Proposition \5.5\ The proof follows naturally from geometric considerations. 
Proof of Lemma \5.3l Letting \t\ < £ + 1/2 and recalling equation (|18p . we see that 



(26) 



log 



sinc(-7rt) 



E 



i2k 



' k(2k- 1) (£+1/2) 



2fe-l 



oo oo I 



1 



j 2fc (2&- l)(£ + l/2) 



2fc-l 



2 A.: 



k 



k=l j=£+l 

Applying Proposition 15.51 to the function f(t) = when k > 1, we obtain 

-fe A 1 1 



< 



r - 



2(^+1/2)^+1 



(2fc - 1)(£ + l/2) 2fc ~ 1 



Equation (I26p becomes 



E 



2(£ + l/2) f-J V£ + l/2 
fe=i 



2fc , /sinc(-7rf) 
< - lor ' 



QiAt) 



E 



< 0. 



4.2k 



'fc(2fc-l) (£ + 1/2) 



2fe-l 



< 0. 



Restated, 
(27) 



£+1/2 



2(£ + l/2) 1 



FPT72 



+ E 



fc=7V+ 



^(2*-!) (£ + 1/2) 



2fc-l 



< - log 

oo 

s E 

fc=AT+l 



sinc(7rt) 
Qi,K*) 
1 



iV 

E 



,2A- 



J fc(2Jfe- 1) (£+1/2) 



2fc-l 



/-2A: 



fc(2& - 1) (£ + 1/2) 



2fc-l " 



Exponentiating, 



(28) 



i ( 1+1/2) 

'2(1+1/2) j 



/ P. 



1 



,2A- 



v^JV 1_ 

k{2k~ 



e _fe=JV+l fc(2fe-l) (i+1/2) 
,2k 



< 



Qi,e(t)e ^ fc=1 fc(2fc - 1) (^+V2) 

sinc(7rt) 



\2fc-l 



Eoo 
. ,. k=N+l k(2k-l) (£+1/2)2*:-! 



□ 
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2JV+1 

Let £ be chosen large enough so that A{£ + 1/2) 2N + 2 < £ + 1/2. If £ is large enough, then for any 

2N+1 

t G E iN ,A, t = c(£ + l/2)2iv+2 for some c G [--4,-4]. For such t, inequality ([28]) implies 

I 1 N + 2 -EJ- 1 C 2fc fiM/n-lfl-TAr) 

g\ 2(^+1/2) W+T l- c 2(^+i/ 2 ))iV+l / e 2^fc=JV+l fc(2fc-l) (,«+ 1 / 2 ) V + / 

_ y-JV 1 t 2k i \ 

< gv(^ M2fc ' 1) (w/2)2fc ' 1 < e E^ +1 ^(W 1 "^) 

sinc(-7rt) 

, 2iV+l , 

If t G (E^a) , then i = + 1/2) 2 ^+ 2 for some c G [--4, A) d . For any such i, we have 

( d N+2 =r- ) ^00 M\%1. ft 1 1 /^(i-iAr) 

(29) g\ 2(£+l/2)^+T l~A2(f+l/2))W+T / e 2^fc=JV+l fc( 2 fc-l) V +/ 

1 11*111 „, , \ 

SINC(vri) 

On one hand, 

/ \ / llcl| 2(JV + 1) -1 \ 

(30) gE^+i^W^) < e (nJfel+°((^/ 2 )^)) 
where the "big O" constant is independent of c G [—-4, A] d . On the other hand, 

(31) g(JV+l)(2JV+l) < e Z.fe=iV+l fe(2fc-l) «+l,^ 



Inequality (|29j) yields 

. d A 2 || C || 2(JV+1) 

(32) ( e 2(^+1/2)^ + 1 i_A2(<+l/2))JV+l _ 1 j e (JV+l)(2JV+l) 

_ \piV 1 ll'llifc . „2(JV+1) 

G>Hp(t)e k=1 k{2k ' 1) " C|l 2(iV+l) 

< y g (AT+l)(2iV+l) 

SINC(vrt) 



dA 2 (^+l) / O, - 

< g(AT+l)(2AT+l) I g \(«+l/2) JV+1 

_ JV+2 

The left most side of inequality (|32|) is of the order 0((£ + 1/2) JV + 1 ), and the right most side of 

inequality (1321) is of the order 0((^ + l/2)~ JV + 1 ), where the "big O" constants are independent of 
c G [— -4,-4]^. The lemma follows readily. □ 

Proof of Lemma \5.4\ Equivalently, we need to show 



|,2(JV+1) 
M 2(iV+l) 

lim sup 

e ^°° ce[-A,A] d 



/ l|C|l 2(iV+l) x 2AT+1 , 

(g(iV+l)(2]V+l) _ lj f(p(t+ l/2)2JV+2) 



0. 
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Suppose the contrary. Let q G [—A, A] d be a value that maximizes the £-th term in the above 
limit. There exists (4)fceN; an d e > such that for all k G N, 

ll r |,2(iV+l) 

/ 11 N 2(jV+l) \ 2JV+1 . 

e < sup e (iv+D(2iv+i) - l)/(c(4 + l/2)aw+5) 

dA 2(iV+l) v 2JV+1 

< [eW+W - lj |/(c 4 (4 + 1/2) 2W+2) | , 

/ 2JV+1 

so that the sequence [f[ce k (£k + 1/2) 2JV+2 )) fcgN is bounded away from 0. This implies there exists 

II 2N+1 .. 2N+1 

5>0 such that \\c ik (£ k + 1/2) ^+2 \\ 2{N+1) < 6 for k G N, that is, ||c 4 || 2 (jv+i) < $(£ k + 1/2) 2iv+2 . 
This forces 



e < sup 

ee[-A,A] d 



" "2(iV+l) \ / 2JV+1 , 

(N+1K2N+1) _ lj/(c(4 + 1/2)5^+2) 



e 

5 2(JV+1) 



< ^ e (« fc + l/2)2^+l(JV+l)(2JV+l) _ 1^ J 

The last term in the above inequality has limit as £ — > 00, which is a contradiction. □ 
Now we can prove Theorem 15.11 



Proof of Theorem \5.1[ If / G PWt^ jj, Theorem 14.21 states that 



Qd,£\b) 

where ^ — > on M d both in L2 and senses. By Lemma 15.31 we have 



(33) sup 

te(E e , NiA ) d 



*<(t)e" 



"1 



,,2(Af+l) 



2^=1 fc(2fc-l) (£+1/2)2*-! _ e (f+l/2)2JV+l(iV + l)(2iV+l) SINC(7Tt) 



<M{£ + l/2)-^m S up (|/(*)|-|&(*)l), 

te(^,jv,A) d 



the right side of which has zero limit. Also, 



(34) 



sup 

i6(^,jv,A) d 



" 2 ' iV + 1 ) \ ^(t) 



< sup 

te(£ f ,]v,A) 

/ dA 2 < iV + 1 ) 

| g(JV+l)(2JV+l) _ I 



/ -^v-i-i; \ Welt) 

( e ( f +i/2)2iv+i (iV+ i)(2iv+i) - 1) n / SINC(Trt) 

^ ' Qd,i(t) 

/ ' "2(JV+1) , 

I e (£+l/2)2JV+l(JV+l)(2JV+l) _ 1 j 



+ 



sup 

t6(E^jv,A) d 
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whose right hand side, by Lemma 15.44 also has zero limit. Combining inequalities (|33p and ([3 
we obtain 



lim 

l— >oo 



' 2^k = l 



„2k 
n 2k 



fc(2fc-l) (£+l/2)2fe- 



Equation (123|) follows by a final application of Theorem 14.2 



SINC(vrt) 



L aD ((E e , N , A ) d ) 



Now we prove equation (122p . Lemma 15.31 and Theorem 14.21 imply 

2fc llfll 2 '^ 1 ' 

Y-N 1 |% m '2(N+l) Hfi(t\ 

Z-fc = l fc(2fc-l) (£+i/2)2fc-l _ p (£+l/2)2JV + l(JV+l)(2JV+l) * V 1 



(35) */(t)e 

< + l/2)-^TT ||/ + ^|| L2(( ^ An , 
the right hand side of which has zero limit. Also, 



QdAt) 



SINC(Trt) 



(36) 



„,," 2 ( JV + 1 ) \ \fj„(t) 

(£+i/2)^+i (JV+1)(2JV+1) _ i) *y siNC(Trt) 

ll t |,2(Af+l) 

/ " "2(iV+l) s 

I g(£+l/2)2JV+l(JV+l)(2iV+l) — 1 J /(t) 

Ht|| 2(JV+1) 
, " "2(JV+1) . 

I g (£+l/2)2JV+l(JV+l)(2JV+l) _ l)^ £ (t) 



^2((S£,iV,A) d ) 



^2((S£,JV,A) d ) 



^2((i?£,iV,A) £i ) 

The second term in the right hand side of inequality ([36]) is bounded from above by 

dA 2(N+l) 



(eWTiTTO-lj|| 6 || i2((Bw)d) , 



which has zero limit. The integrand of the first term in the right hand side of inequality (j36|) (as a 
function over M. d ), converges uniformly to zero by Lemma 15.41 and is bounded from above by 



e(N+1)(2 N+i) _ lj|/(t)| 2 G ), 



so this term has zero limit by the Dominated Convergence Theorem. Combining equations ([35 
and ([36]) yields 



lim 



**(t)e" 



!!aL 



fc(2fc-l) (£+l/2)2fc-l 



CM*) 

Equation (j22p follows by a final application of Theorem 14.2 



SINC(TTt) 



L2((E ttN , A ) d ) 



0. 



□ 



The optimal growth for any (-E^jvfc)^ such that Theorem 15.11 holds is not known, but an upper 
bound for the rate is established in the Appendix. 
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6. Examples of uniformly invertible exponential Riesz Bases 

Theorems 14.21 and 15.11 both require uniform invertibility of (fn) n eZ d - Fortunately, there are signifi- 
cant classes of exponential Riesz bases (f n ) n ^z d which have this property. Consider the exponential 
Riesz bases described given by Theorem 16.11 (Corollary 6.1 in [JJ) and Theorem 16.21 (Theorem 1.3 
in [E ;. 

Theorem 6.1. Let (tk)k&z d c ^ d such that 

r ln(2) 
sup \\n - tnWoo = L < — , 

n&L d ncl 

then the sequence (fk)kei, d defined by fk{x) = / 27r W2 e i ^ x ' tk ' is a Riesz basis for L2Q— 7r, 7r] d ). 

Theorem 6.2. For d > 1, define 

_ , . / . sin7rx\ rf /sin7rx 
Dd{x) := 1 — cos irx + sm ttx H — 

\ TTX ' \ TTX 

and let Xd be the unique number such that < x& < 1/4 and Dd(xd) = 1. Let {tk)kez d ^ ^ d such 
that 

SUp \\n - tnWoo < X d , 

nei d 

then the sequence (fk)kez d defined by fk(x) = ( 27r ^/2 e l( - x ' tk ^ is a Riesz basis for L 2 ([— vr, Tr]^). 

It is worth noting that when d = 1, Theorem 16.21 reduces to the classical Kadec's 1/4 Theorem, 
first proven in [3]. A proof of Kadec's 1/4 Theorem can also be found in |13} Theorem 14, page 
36]. The proofs of the above theorems show that the map Le n = f n satisfies || I — L\\ = 5 < 1, 
from which we see that L is invertible. Uniform invertibility is readily verified, or can be seen as a 
consequence of the following more general proposition. 

Proposition 6.3. Let L : H — > H be a uniformly invertible operator with respect to (Pi)g^fq, where 

limsuplKPfLP^ 1 !! = M < oo. 

I— >oo 

Lf A is an operator such that 

(37) \\L-A\\ = £ 

for some 7 < 1, then there exists N £ N such that A is uniformly invertible with respect to (Pi)i>n. 

Proof of Proposition \6.3i Using uniform invertiblity of L, and noting that (L^) _1 — (L*) _1 = 
(L^) _1 (L* — L^)(L*) _1 and lim^oo L*^x = L*x for all x £ H , we see that lim^_ 5 . 00 (L^) _1 2; = (L*)~ 1 x 
for all x 6 H. Equation ([6]) implies that lim^ 00 (L|) _1 x — (PeL* P()~ 1 x = for all i£ff. Together 
we have 

lim (P t L*Ptf x x = (L*)~ 1 x, xeH. 

£— »oc 



General principles imply 



< liminf UP.LPtYH < M. 



I— >oo 
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This yields \\L — A\\ < uj^-iu ; so that A is invertible by usual Neumann series manipulation. 

Equation (|37p yields that \\PgLPt — PgAPgW < jj, where the norm is now the standard matrix norm 
on the set of matrices of dimension dim(ranP^). This implies 

\\Pi - {PiLPty^PtAPt)]] < ^IKP.LP,)- 1 !!, £ e N. 
Choose N large enough so that || [PgLP^ 1 1| < when £ > N. This yields 

\\Pi - (PiLPt^iPtAPt)]] < 2±l t £ > at. 
Standard manipulation shows that P^APi is invertible for £ > N, and 

sup UPiAP^W < Z +1 , M. 

£>N 7(1-7) 

□ 

Note that in the previous proof, if M is redefined to be sup^ gN || (i-^-Li-^) - 1 1| , then A is uniformly 
invertible with respect to (Pi)i^. 

The following proposition shows that compact perturbations (of arbitrary norm), of a uniformly 
invertible operator also gives a uniformly invertible operator. 

Proposition 6.4. Let L : H — > H be uniformly invertible with respect to (Pi)?^. If A : H — )• H is 

compact such that L = L + A is an onto isomorphism, then there exists N such that L is uniformly 
invertible with respect to (Pe)i>N- 

Proof of Proposition \6.4\ From the definition of Lg, we have I = (I — Pe)Lj 1 + LPpLj 1 , so that 
This implies 



L-\P £ - I)LJ X = P t L- x - L-\ 



(38) m~ x P t - (L*)- 1 = {L\)-\P t - I){L*)-\ 

As £ — > 00, the right hand side of equation (I38p has limit pointwise. Combined with the com- 
pactness of A, we obtain 

lhn(L})- 1 P e A = (L*)- 1 A 

£—>oo 

where limit is in the operator norm topology. This yields 

(39) lim / + APiLj 1 = I + AL" 1 = (L + A)L~ 1 , 

where the limit is also in the operator norm topology. The right had side of equation (j39|) is 
invertible, so there exists N such that £ > N implies (/ + APgLj 1 ) -1 exists, and 

(40) sup ||(/ + AP^ 1 )- 1 !! < 00. 

e>N 

From the definition of L^, we obtain 

L e = L e + AP e = (I + AP e Lf)L e . 

When £ > N, we have 

LJ 1 = LJ 1 (I + AP e LJ 1 )-\ 
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and equation (|40p implies 



(41) 



sup \\L 

t>N 



-li 

I 



< sup \\L 

£>N 



-li 



sup 

e>N 



(I + AP.L,- 1 )- 1 ! 



< oo. 



This completes the proof. 
The following lemma holds: 
Lemma 6.5. Choose (tk)keN C 



□ 



k=l 



such that (h k ) ke n ■= ( (27r ) d / 2 e' ((,),ffc) ) fcgN satisfies 

1/2 



< 



L 2 [-7T,7r] 



fe=l 



//(r fe ) fc6 N C and (/ fc ) fc6N := (^-^(('W) 



for all (c fc )Li C C. 
£/ien /or a// r, s > 1 and any finite sequence 



(o>k)k, we have 



k=r 



< Be 



L*2 [ — TT,7v] d 



j ( sup 

7rrf I , , ^ 

3 \r<fc< 



1 



k=r 



This lemma is a slight generalization of Lemma 5.3, found in pQ using simple estimates. Lemma 
6.51 is proven similarly. A consequence of Lemma 16.51 is the following corollary. 



Corollary 6.6. Let {tk)k& C M rf such that (/fc)fceN (the usual exponential sequence defined in 
terms of (t n ) n ) is a Riesz basis for 1*2 ([~ tt] )■ (ffc)fceN C M d , and define (gk)k£N by gk{%) = 
L e'^'^l Lei (6jfc)fce2 ^ e an orthonormal basis for L 2 ([— vr, 7r] d ). 7/ 

(27TJ 



lim 



7~fc 1 1 oo — > 

i/ien i/ie operator defined by Ofc i— >■ — 5% *s compact. 

The proof of Corollary 16.61 is similar to the proof of Corollary 5.5 in [1], so it is omitted. 

Corollary 6.7. Let (ifc)fceN C Let (//OfceN be a Riesz basis for L2Q— 7r , ir] d ) which is uniformly 
invertible with respect to the projections (Pe)een defined at the beginning of section 4- If ( T fc)fceN C 
M d , and (gk)kef$ are as * n Corollary W. 6\) . and additionally, (gk)k&N is a Riesz basis for L2([— 7r,7r] d ), 
then (gk)k£N is uniformly invertible with respect to a sequence of projections (Pe)e>N for some 
N > 0. 

Proof of Corollary \6. 1\ Apply Corollaries 16.41 and 16.61 □ 

This corollary relates to Theorem 14.21 in the following way. Let (gk)keN is as in the preceding 
theorem. Usage of Corollary 16.41 in the proof of Corollary 16.71 does not ensure that low order 
polynomial interpolants will exist; however, they will existence for sufficiently large I. Simple 
examples show that in Corollary 16.71 the additional assumption that (gk)keN is a Riesz basis for 
p2([— 7T, 7T ] d ) cannot be dropped when d > 2. Example: The standard exponential orthonormal 
basis {e n )n£_i d ls °f course uniformly invertible, but the set 

I J((-W,V2fl,- ,o))\ u ( \ 
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is not a Riesz basis. However, this condition can be dropped when d = 1. This follows from Corol- 
lary [6]6] and the following theorem. 



Theorem 6.8. Let (i n )nez C M be a sequence such that (f n )nez (defined as before) is a Riesz 
basis for L 2 [— vr, n]. If (r n ) ne 2 C R is a sequence of distinct points such that 

lim \t n - T n \ = 0, 

\n\— »oo 

then (g n )n& (defined as before) is a Riesz basis for L 2 [— 7r,7r]. 

The proof of Theorem 16.81 relies on Lemma 16.91 below, which originally appears as Lemma 3.1 in 
[7J. The proof of Lemma 16.91 found in jJJ itself relies on a citation, so for the sake of completeness 
it is presented here with a self-contained proof. 

Lemma 6.9. Let (f n )nei< be an exponential Riesz basis for L 2 [— 7r, 7r]. If (g n )\ n \<e * s a sequence of 
complex exponentials such that 

(g n )\n\<e U (fn)\n\>£ 

consists of distinct functions, then this set is a Riesz basis for L 2 [— 7r, 7r]. 

Proof of Lemma \6.9\ If we can prove the case when I = 0, the general result follows inductively. 
Let /„(•) = ie'*' 1 ") for and g (-) = -^e 1 ^ where r G R and r / t n for n ^ 0. Let 

(e n ) n £z be an orthonormal basis for L 2 [— vr,7r]. Lemma f3.5l shows that (go) U (/ n )n^o is a Riesz 
basis if and only if the map defined by 

eo ^ 5o, e k i-> f k , for k / 

is one to one. This is readily seen to be equivalent to (<70)/i) 7^ 0, or by passing to the Fourier 
transform, to Gq(to) ^ 0, (recall that Go = J-(f^)). If we can show that the only zeros of Go in R 
are (i n )n^o> we are done. 

Suppose there exists A G R, A ^ (in)n^o- Such that Go(A) = with multiplicity m. Define the 
entire function 

*<«> - fF^>- 

Note that H\r G L 2 (R), and H is of exponential type tt, so H £ PW\- n ^-\ by the Theorem 12.41 
The expansion 

H(t) = Y, H (tn)G n (t), 

combined with H(t n ) = 5 n fl, shows that H(t) = Go(t) for all t G R, an immediate contradiction. 
We conclude that G (A) ^ 0. □ 

Proof of Theorem \6.8[ Define Le n = f n and Le n = g n . By Corollary 16.61 L is bounded linear and 
L = L + A for some compact operator A. Define the operator 

fn, \n\ < i 
9n \n\> £ 



Rewritten, we have 



R £ = LP £ + (L + A) (I -P e ) = L + A(I - P ( 



e 
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Compactness of A implies that lim^oo = L in the operator norm topology. We conclude that 
Rl Q is an onto isomorphism for some £q sufhcently large; that is, the set 

( 42 ) (fn)\n\<£o U (9n)\n\>£ 

is a Riesz basis for L%\— tt, 7r]. If we apply Lemma [6.91 by replacing (f n )\ n \<e with (g n )\ n \<e m 
expression ()42|) . we have that (g n )nez is a Riesz basis for L2[— vr,7r]. 



□ 



7. Appendix: Comments regarding the optimality of Theorem 15.11 



In the statement of Theorem 15. 1\ it is not apparent whether or not (En N,k)e can be replaced with a 
more rapidly growing sequence of intervals; however, Proposition 17. II shows that if f(t) = SINC(7rt), 
equations (f23|) and ([22]) can hold for a sequence of intervals [E^)i which grow faster than (E£^,A)e- 
Propositions 17.61 and 17.81 show that growth bounds of the intervals in Proposition 17.11 are optimal 
for the conclusion of said proposition to hold. Thus, the bounds in Proposition 17.11 provide upper 
bounds for the growth of any sequence (E^^^t such that either equation (|23l) or equation (|2~2j) 
hold for general multivariate bandlimited functions. 



Proposition 7.1. Define 
CfN = 



D 



LN 



-(2N + l) 2 (£ + l/2) 2N+l \og{l + 1/2) 
\{2N + l) 2 (£ + l/2) 2N+l log(£ + 1/2) 



i 

2(]V+1) 



1 

2(JV+1) 



and 



then 
(43) 

and 
(44) 



lim 

l— >oo 



lim 



SINC(TTt) - Ismcn(.)At) 



SINC(vrt) - IswcnuM 



L ao (l-D etN ,Dt, N ] d ) 

The proof of equation ()43|) requires the following two propositions. 
Proposition 7.2. 



0. 



(45) 



lim 

£— >oo 



llt|| 2(JV+1) 

f 11 N 2(JV+1) 

g (£+l/2) 2JV + 1 (iV+l)(2iV+l) 



1 SINC(vrt) 



0. 



Proof of Proposition 7.2. Let t = aCe t N where a G [—1, 1] . Noting that 

1 2jv+1 llall 2(JV+1) 
y£+ 4(JV+1) 11 " 2 ( JV+1 ) 

the quantity in equation (|4"5j) becomes 



H,,|2(iV+l) 
l|c|l 2(iV+l) 
g (£+l/2) 2A ' + 1 (]V+l)(2]V+l) 



< 



[— Ce,N ,Ce y N] d 
1 



+ 



1 



4(]V+1) II ll2(JV+l) 



2JV+1 



a 



d/2 
LN 



-1,1]° 



4 



2JV+1 || ||2(JV+1) 
4(AT+1) II H2(Af+l) 



1 SINC(TTt) 



2 v 1/2 

da I 



dt 



1/2 
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< 



2-/2^+1) ™ +2^/2 



(|(2iV + 1)2) 4(iv+Ty ( log (£ + 1/2)) 4(iv+D + i/2) d 3(w+ij 



The last term in the above inequality has limit as i — > oo. This proves the proposition. □ 
Proposition 7.3. 



lim 



,2k llfll 2 < JV + 1 ' 

_ r H 1 Wife ' "2(iV+l) 

Qde(t)e ^ k=1 (/+l/2)^-l_ e («+l/2)^+l(JV+l)(2JV+l) SINC(TTt) 



= 0. 

ia([-C.jvA,Jv] d ) 



Proo/ o/ Proposition\7.3\ If i G R d and p||oo < £ + 1/2, then equation ([28]) implies 

/ \ 2 



/ y^oo 1 11*111 \ * 

(46) ( e ^k=N+l k{2k-X) (£+1/2)2^-1 ) TT 



v» i \ " \ 2^+1/2) / \ 

z^fe^iv+i fc( 2fc-i) (f + i/ 2 )2fc-i i TT e v ^Ki+tjz) 



1 V 1+1/2, 



fc=i 

i2fc 



O, Jf)p ^fc=l fc(2fc-l) (/+l/2)2fe-l 1 W|fc 

< < 2^ fc=JV+ i fc (2fc_i) (f + i/ 2 )2fc-l 

SINC(Tri) 

Let i E [— C^jv, Q,Ar] d where i = aC^, a G [—1,1]. Consider the right hand side of inequality 
I6l) for such i. 



^oc i Ml / 1 \ 2JV+1 ||a|| 2(JV+1) r£+l/2lofll * || 2(JV+2)> 1 

(47) e E fc=J v+i I(2Fr7y ( , +1/2) 2fc-l < ^ + i .^4(iV+i)ll Q ll2(iV+i) e ^+l/^W^|l+T72|l2(iv + 2); 

-, 2iV+l I, m2(JV+1) 1 N+2 , , 

< (l + I^5(iV+T)!l a ll2(iV + i) e M(^+l/2) iV + i ( i og( ^ + i /2 ))7V+T || Q |j 

for some constant M. Noting that 

_ ||Q|||(|(2jV + l) 2 )^T(log(£ + l/2))^I 

(^+V2) 3_ (/+i/ 2 )im 



we can bound the left hand side of inequality (|46h from below as follows: 



(- 



i 

«2l ,„ r .,,^ m )^TT 



Hallo I log(£+l/2) 



TV+S I , 1 ■. 2Ar+1 HaH 2(iV+1) 

(48) e V (^i/2)*ff / ( i+ l^iW+IjI^CJv+i) 



/ ^oo 1 W|| \ ^ I 2(£+l/2) / t \ 2 I 

< I e 2 - fc=JV + 1 fc(2fc-i) (£+i/ 2 )2fc-i I TT e 1_ lw72/ 7 

\ 7 ; — i 



k=l 

where m > is chosen independently of I. Relations (|46p through (|48p imply 



j 

~~ ~ 1 1 2 



[e K <^/2)^t / 4(JV+1) " ' l2(JV+1) |SINC(vrt)| 
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< 



at 1I + I| 2(JV+1) 
_v-JV 1 Plljjj |f|l 2(iV+l) 

Qde(t)e ^ k=1 k{2k - 1) (^+i/2) 2fc - 1 - e (^+1/2)^+1 (JV+ i)(2iv+i) SINC(Trt) 



< e 



M- 



(log^+1/2)) 



lVTT|| a |,2(Af+2) 
11 "2(iV+2) 
1 



(e+i/2) TO 



2AT4-1 ,,2(JV+1) 



4(]V+1) ll""2(iV+l) 



|SINC(vrt)|. 



Further simplification implies (for appropriate constants C, C", and C") that 



ii 2 *' 



imi 2(JV+1) 

11 "2(iV+l) 



Qd,i{t)e ^ h=1 k ^ k - 1) (t+VV™- 1 _ e (f+i/2)2«+i(jv+i)(2iv+i) siNC(Tri) 



< C 



JV+2 

(log(£ + l/2))^+i 



(£ + 1/2) 



JV+l 



C 



< c"- 



JV+2 

,,(log(£ + l/2))*TT 



+ 



«*5 



2N+1 i n 2(JV+l) 



L 2 (l-C e>N ,C e , N ] d ) 

2 \ 1/2 



4(JV+1) H""2(JV+1) 



x 2JV+1 || i|2(jV+l) 
]\ 4(JV+1) ""H 2 (JV+1) 



a|||SING(7rt) 

2 



a||!SINC(7rt) 



dt 



1/2 



-da 



(£ + 1/2)^+1 \J[-i,i]* (\og(£+ 1/2))2W+TJ^ + i)2Hv+ry« 



iV+2 

(log(£ + l/2))^+T 



(£ + l/2)5v+T(log(^+ 1/2))^+^ 
after the change in variable t = ctCi^ and simple estimates. This proves the proposition. 

Proof of equation ft43\ ). Equation (]43p follows immediately from Propositions 17.2 1 and 17.31 

The proof of equation (I44D requires the following two propositions. 

Proposition 7.4. 



□ 
□ 



(49) 



lim 

i— >oo 



,2(iV+l) 



^ e (£+i/2)2jv+i ( jv+i)(2jv+i) _ i"jsinc(7ri) 



0. 



Proof of Proposition \ l.J\ Let t G [— D^jsii D^jy], then t = aDi^ f° r a £ [ — 1>1]- Simplification 
shows that equation (j4"9|) holds if 



(50) 

Note that for large £, 



lim sup 

e->ac ae [ ,i] 



„2(JV+1) 2JV+1 

(£+1/2) 2 < Ar + 1 ) - 1 



2N+1 



a 



( log(£ +1/2)) 2 ( Ar +D (£ + 1 /2) ^W+T) 



0. 



ri ,2(JV+l) 2JV+1 
(£+1/2) 2(JV+1) _ i 



(51) sup 

ae[i/2,i] a (log(£ + 1/2)) W+T)(^ + i/2)5jwtfy 

Let < a < 1/2. The Mean Value Theorem implies 



< 



log(£ + 1/2)) 2(~+D 



(52) 



„2(JV+1) 2iV+l 
(£+1/2) 2(iV+T) _ i 



a 



< (2JV + !)(£ + 1/2)° 3P^a 2JV+1 log(£+l/2). 
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This yields 



sup 

<*e[o,i/2] 



n 2(N+l) 2JV+1 
(£+1/2) 2(AT+1) _ 1 



2JV+1 



< M- 



log(£ + 1/2)) 



(1 + 1/2) 



^±^(1^(1/2)2(^+1)) 



a(log{£+ l/2)) 2 (^+D (£ + 1/2)^+^ 

for some constant M. Combined with inequality (|5ip . we have equation (|50p . which proves the 
proposition. □ 



Proposition 7.5. 



lim 

t— >oo 



t 2(N+l) 



_ y^JV 1 t ZK 

Qlt(t)e ft=1 fe(2fe - 1) (f+i/a) 2 *- 1 - e («+i/2)2iV+l (A r +1)(2A r +1) ging^) 



0. 



Proof of Proposition \ 7. 5\ Let i G [— C^jv, C^,jv] where i = aC^jy, a 6 [—1,1]. Proceeding in the 
same manner as in the proof of Proposition 17.31 we see (for appropriate constants C and C") that 



V JV 1 t 2(iV+l) 

L-*k=\ fc(2fc-l) (<! + l/2)2fc-l _ e (i+l/2)2JV+l (Ar+ i)(2Ar + l) 3^(7^) 



< 



< 



„,2(JV+1) 2JV+1 AT+2 

1/2)° ^+T)q2( log (£ + l/2))^+T|sin(7rt)| 

1 1 2iV+l 

a(£ + l/2)^+r(log(£ +1/2)) (£ + 1/2) ajw+lj 

2iV+3 

C"(log(£ + 1/2)) 2 (^+D 



(£ + 1/2)^+1 
This proves the proposition. 



□ 



Proof of equation {4$ - The previous two propositions prove equation when d = 1. The multidi- 
mensional case follows inductively. □ 

Proposition 7.6. Let N > 0. // (M^jv)^ is a sequence of positive numbers such that |^3| ) ZioWs 
w/ien (Ci t N)e is replaced by (M£^)i, then 



(53) 



limsup^+ — < 1. 



N 



The proof of Proposition 17.61 requires the following simple estimate: 
Proposition 7.7. Let a > 1/2, e > 0, < u < 1, then 

e a 



(l+e)a gin 2 ^ x 



,1+a; 



dx > 



2a w (l + e) w 2(a-l/2) 1 +" 



Proof of Proposition \ 7. 7\ Let b = (1 + e)a. We have 



■ 2 

Sin 7TX 



dx + 



2 

COS TTX 



„1+UJ 



H 1 

u}\a u 



and 



b i 
COS 7TX 



n.l+Ul 



-dx 



b - 1 ' 2 sin 2 7rx 



<ix < 



0-1/2 (X+1/2) 1+W 1-1/2 



^V2 8^2^ 



-dx. 
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This yields 



so that 



Noting that 



b ■ i 
sin 7TX 



dx 



b ■ 2 

sin TTX 



6-1/2 



sin 7tx 1/1 1 

r- OX > — — 

a- 1/2 



6 sin 2 ttx 1/1 1 
ax > — I 

x 1+ul 2ui \ a" 



1 



1 / 1 



1 



2ta Va^ b u 
proves the proposition. 



2(a - 1/2)!+^' 
(l + e) w -l 



2wa a; (l + e) 



> 



2a w (l + e) 



□ 



Proof of Proposition 7.6. Fix ./V > 0, and define c = 2N+4 where < o is small enough so 

that c < 1/2. Define 



and 



^ = ( c (iV + 1)(27V + l)log(£ + 1/2))^+^ + 1/2) ^W+iJ 



e^^ + 1/2) 1 - 2 ^. 

Note that lim^^ e £ = 0. Let i G L4^ + 1/2), (1 + e^)^(£ + 1/2)], then t = a{i + 1/2) for some 
a G L4^, (1 + t()Ag\. For large £, note that inequality ([28]) implies 



, , 1 «+i/2)a a ^ +1 ) Isin7ra(l + 1/2)1 
(54) — e (iv+i)(2iv+i) J i — ! — !—L± < 

1 j 2vr a(£+l/2) ~ 



Ql,l(t)e 



v^JV 1 



fe(2fc-l) (£+1/2)2*- 



Moving to the multivariate case, if t G \A t {£ + 1/2), (1 + e e )A e (£ + l/2)] d , then t = a(£ + 1/2) for 
some a G [Al, (1 + e,)A,] d . This yields 



n 



1 |sin7rai(^ + 1/2)| 
JM (^(^ + 1/2))^ 

For sufficiently large £, we can conclude 

1 



< 



Qdl(t)e 



Hi 



fc(2fc-l) (£+i/ 2 )2fc-l 



(1+^)^+1/2)^2^ id 



9vr 2 Af 7A,(€+l/2) 



< 



„2-2c 



■dx 



Q<Le(t)e 



fc(2fc-l) (£+i/ 2 )2fc-l 



'[ J 4 £ (£+l/2),(l+e < )AH<?+l/2)] d 

Applying Proposition 17.71 for a = + 1/2), e = e^, and a; = 1 — 2c, and using the definition of 
eg, we obtain 



l r 



1 



1 



9tt 2 L2(l + e^)!- 2c 2^ 2c (^(£ + 1/2) - l) 2 ~ 2c 



< 



[A t (t+l/2),(l+ei)Ai(i+l/2)]* 



Qd,e(t)e 



fc(2fc-l) (£+l/2)2fe- 
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The first term in the brackets in the previous equation has limit 1/2, while the second term has 
limit 0. We conclude there exists a constant (3 > such that 



(55) 



Z^fc = l 



dt, £>0. 



l[A e (e+l/2),(l+e e )A e (e+l/2)] d 

If Mt jq > (£ + 1/2)(1 + eg)Ag for infinitely many £, there exists a subsequence (£k)keN such that 
(in particular), 



lim 



SmC(7rt) - Q dA (t)e 



Z^fc=i fc(a*-i) (i fe+ i/2)2fc-i 



ia([^< fc (4+l/2)),A^ (4+l/2)(l+e^ 
This contradicts inequality ([55]) . This yields that for sufficiently large 
M £iJV < (£ + l/2)(l + e e )A e 
'2N + 1 



■ {). 



1 + Q 



4iV + 4 



+ (5/2 ) (2V + 1)(2N + !)(£ + 1/2) 2N+1 log{£ + 1/2) 



2(iV+l) 



2(JV+1) 



Note that since eg — > 0, for large £, the quantity (1 + ( 4^+4 + ^/^ 

bounded away from) the quantity (f^j + 5 
£n,s such that 



1 



2(JV+1) 



is less than, (and 



. We conclude that for any 5 > 0, there exists 
2iV + 1 



sup 



2JV+1 ^ ( 
1/2)2(JV+1) V 



4iV + 4 



+ 8 



2(iV+l) 



□ 



*>4m ((AT + 1)(2JV + 1) log(£ + 1/2)) 2 ( N + 1 ) (£ + 1/2) 
Proposition 17.61 follows. 

Proposition 7.8. Let N > 0. If {Mg^)g is a sequence of positive numbers such that equation fi44\ ) 
holds when (Di^e is replaced by (Mi^i, then 



(56) 



V M ^ N <: 1 
limsup — < 1. 



e^oo Din 

The proof of Proposition 17.81 requires the following fact: 

Proposition 7.9. Let < e < 1. If A > 0, there exists t G [A, A + e] such that |sin(7ri)| > 
I sin(vre/2)|. 



Proof of Proposition [7T9[ The proof is clear from geometric considerations. 

Proof of Proposition \Z8\ Let N > 0. Choose 5 > such that c := |^±| + 5/2 < 1. Define 



and 



□ 



A e = (c{N+ 1)(2N+ l)log(£+ 1/2)) 2(^+D (^ + 1/2)2(^+1) 



e e = A e (£ + l/2) 1 - c . 

Note that lim^oo eg = 0. Let t S [Ag(£ + l/2), Ag{£+l/2) + eg}. Proceeding as before, for sufficiently 
large £, we have 

t 2(JV+l) 



_ A (e+l/2)™+l (N+1 )(2N+l)) 

2ir 



sin(7rt)| 
t 



< 



Ql,i(t)e 



EN 
k=l 



fc(2fc-l) (f+l/2)2fc-l 
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Now for all t G [A e (£ + l/2),A e (£ + 1/2) + e e ], 



1 (l + l/2) c ,„ : „^, ^ ^ -£JU 



sin(7rt)| < Qv(t)e ^ ft=1 fc(2fc - x) c+i/2)«-i 



2tt + 1/2) +e/ 
In the multivariate case, if t G [At(£ + 1/2), Ai(l + 1/2) + e^] d , we obtain 

1 (C + \l2Y d -A- -v^ i Mil 



(2Tr) d (At(e + 1/2) + e e ) d f \ 

For large £, an application of Proposition 17.91 yields 

1 |sin(W2)| rf „ _i — 

< Qd,i{t)e k=1 fc(2fc " 1) (^+i/2) 2fc - 



(3vr) d ^(£ + l/2)( 1 - c ) d 



ioo([^(^+l/2),^(^+V2)+Q] d ) 



By the definition of q, the right hand side of the above equation tends to a positive constant. The 
remainder of the proof is almost identical to that of Proposition 17.61 □ 

The following is trivially deduced from Propositions 17.61 and 17.81 Fix N > 0. If (E^at^ is a sequence 
of intervals such that either equation (|23p or equation (|22|) holds for all / G PW^_ n ^d, then 



max \\x\\oo = °l max | x||, 
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